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Q^' The space of local operators in massive deformations of conformal field theories 

. is analysed. For several model systems it is shown that one can define chiral sectors 



in the theory, such that the chiral field content is in a one-to-one correspondence 
with that of the underlying conformal field theory. The full space of operators con- 
sists of the descendent spaces of all scalar fields. If the theory contains asymptotic 
states which satisfy generalised statistics, the form factor equations admit in ad- 
dition also solutions corresponding to the descendent spaces of the para-fermionic 
operators of the same spin as the asymptotic states. The derivation of these results 
uses g-sum expressions for the characters and g-difference equations used in proving 
Rogers- Ramanuj an type identities. 
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1 Introduction 



During the last several years the theory of completely integrable models in (1 + 1) di- 
mensions has experienced a rapid development. Especially the study of conformal field 
theories (CFT), which resulted in the solution of several classes of such theories, has 
allowed a deep understanding of the critical points describing second order phase transi- 
tions. Unfortunately, less is known for off-critical integrable field theories. Their solution 
should give a description of the scaling region around critical points and allow a deeper 
understanding of standard techniques used in quantum field theory, such as the renor- 
malisation group or perturbation theory. 

In this article we will study so-called perturbed conformal field theories, i.e. models 
which can be viewed as perturbations of a conformal field theory by a relevant operator. 
The perturbation is chosen in a way that the theory remains integrable also away from 
the critical point. This description proves useful, since many properties of the conformal 
theory can be preserved or 'deformed' off-criticality. As, for example, the quantum group 
structure, which becomes an affine one in the perturbation process, or the conserved 
charges. 

The ultimate goal in the investigation of these systems is to solve them, that is 
to determine the algebra of local operators and their correlation functions. In general 
conformal field theories perturbed by the operators 0i^2, 02,i and 0i^3 

H = HcFT + d^x(j){x) (1.1) 

are driven into a massive regime. It has been shown by A. Zamolodchikov |jl|, that such 
theories can be described by the exact scattering matrix. A lot of progress has been made 
in determining in general the on-shell structure and understanding its algebraic origin. 

The understanding of the off-shell structure of these massive integrable theories is less 
complete. Knowing the exact S-matrix of a theory one can calculate the matrix-elements 
of local operators, called form factors. They are determined from the knowledge of the 
on-shell structure of the theory, through the form factor equations. This is a remarkable 
property, since the on-shell structure of the theory does not contain any information 
about local fields. In fact it is described in terms of asymptotic states. 
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The form factor equations do not refer to a specific operator. This on the one hand 
makes the identification of operators a difficult task. Apart from the fundamental fields of 
a theory (as the energy-momentum tensor or conserved currents) which have additional 
properties, it is difficult to identify a specific operator from the knowledge of its form 
factors. This is true especially in this case, where we are interested in the comparison 
with the conformal field theory in the ultraviolet limit. In order to achieve such an iden- 
tification, one would need to calculate correlation functions, which involves the summing 
up of the form factor series, a problem which has not been solved for interacting models. 

The above mentioned arbitrariness in the form factor equations allows the determi- 
nation of the space of local operators of the theory. Since the equations must be valid 
for all local fields of the theory, one can specify the field content by determining the 
independent solutions of the form factor equations. In this article we will use precisely 
this fact, and find the space of local operators for some simple massive theories. They 
have the distinguishing property that they contain only scalar asymptotic states, and can 
be described as perturbed minimal conformal field theories. 

Beforehand, it is not clear what kind of structure one should expect to find in the 
perturbed conformal theories. Since the perturbation theory defined by ( |1 . 1|) is super- 
renormalizable, only a finite number of counter terms is involved in the renormalisation 
procedure, which should leave unchanged the structure of the operator algebra of the 
theory On the other hand the most salient feature of the critical theory is its Virasoro 
algebra symmetry, which is broken in the massive theory and one expects drastic changes 
in the off-critical theory. In the critical theory left and right chiral sectors (corresponding 
to the components T and T of the energy momentum tensor) are independent which 
certainly is not true in the deformed theory. Also the structure of the conserved charges 
changes in the course of the perturbation. 

The result of this study will be that the structure of the massive theory resembles much 
the one of the conformal theory. We will define a grading in the space of operators which 
allows a formal definition of chiral components. The so-defined chiral structure of the 
deformed theory becomes isomorphic to the conformal one. In CFT on the real line there 
are no constraints on which combinations of the two chiral dimensions of the operators 
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can enter the partition function, apart that they should form a closed operator algebra. 
In the massive model we have an additional locality requirement, deriving from the form 
factor equations. As will be shown in the following sections, only specific combinations 
of chiral dimensions are present in the massive theory. 

1.1 The Form Factor Bootstrap Approach 

Let us review some general properties of the form factors 0, ^ ^, ^. We parametrise the 
momenta in terms of the masses of the particles rrii and the rapidity variables jSi as 

p° = rrii cosh (3^ , p] = rrii sinh [3i . 

Form factors are defined as matrix elements of local operators O between the vacuum 
and the set of asymptotic states, 

•^e?,..,...,.„(A,/?2,...,/3n) = (0|O(0,0)|Z,,(/3i),Z,,(/32),...,^.„(/3n))in . (1.2) 

Properties deriving from basic principles of quantum field theory, such as crossing and 
unitarity, and the properties of the Fadeev-Zamolodchikov operators Z^{(3) allow the 
construction of these matrix elements. For the case of scalar particles these functional 
equations known as the form factor axioms are: 

1) Watson's equations 

•^e*^...e,e,+i...£„(A' • • • ; A, ■ ■ ■ (^n) = 

Sei€i+i{Pi - Pi+l) ^^^...ei+iei...e„iPl^ ■ ■ ■ ^ Pi+l^ Pi^ ■ ■ ■ Pn) , (1-3) 

+ 2vr2,/32, ...,Pn) = e^-C^^+^^+E -^.^ J . . . A) , (1.4) 

where s^^ denotes the spin of the particle e„ and denotes the mutual locality index 
between the field O and the field creating the asymptotic state Z^. Similarly Weie„ denotes 
the generalised statistics, which the operators creating the particles i and n satisfy. It is 
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determined by lim/3_too Sni{f3)- This phase factor is necessary since for some theories we 
will consider, the particles will not become free particles for (3 ±00, but will satisfy 

Therefore whenever one continues the form factors analytically, one needs to take these 
commutation properties into account. For theories containing degenerate particles this 
statistics is governed by the i?-matrix defining the ultraviolet limit of the scattering 
theory Alternatively these phases can be taken in account by assigning non-trivial 
charge conjugation to the operators creating the asymptotic states. 

2) Relativistic invariance 

J^Z.-eM + A, . . . ,/5n + A) = e^'^J^Z.-eM^ ■ ■ ■ , f^n) • (1.5) 

3) Asymptotic behaviour 

J-,,,...,,„(/3i + A,...,A + A,A+i,...,/3n) = 0(e^"l^l) for |A| ~ oo , (1.6) 

and S = max(S'^) is a finite constant. 

4) Equation of kinematical poles 

-t hm (/?' - P)J'Z,...JP' + P,Pi,---, Pn) = 

(l - e2-(-^+S-..) n 5, J/3 - A)) ^^....„(/3i, . . . , /3n) . (1.7) 

5) Bound state equation: Let particles Aj, Aj form a bound state A^, the corresponding 
two-particle scattering amplitude exhibits a pole with the residue 

-I hm {(3-ivZ.^S,.em = ^%^ ; (1-8) 

then 

-i hm (/?' - /3)^„„^,^...,JA, ...,/?' + zu^^,^, /5 - ^n^;,^, . . . , = 

= r^f.,^L......„(A,...,/?,...,/3n-i) . (1.9) 
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In it was shown that the operators defined by the form factors (|1.2| ) satisfy proper 
locahty relations provided that they satisfy the properties ( p..3|) -( pr6D and the residue 
equations ( |1.7| ) and (|1.9D. 

A well established solution method uses the parametrisation of the form factors in 
terms of the minimal two particle form factor and the explicit prescription of the pole 
structure. As a first step one determines the minimal two particle form factor, which 
satisfies Watson's equations for n = 2, 

J'TiPl2) = Sat,iPl2)J'ri-Pl2) , ^r"(^^ - M = (^Vr + A2) , (1.10) 

where Pu = /3i — /?2- In order to specify an unique solution to these equations one imposes 
further restrictions, namely that _F™™ is analytic in < Im/5 < vr and has no zeros in 
this range. The n-particle form factor can then be parametrised as 

^ £l,...,e„(A, • • • , Pn) = 



i<:j {Xi + Xj) '^"^ Yla (smh ^{Pij - la) smh ^[jJij + la^ 

(1.11) 

The product over a runs over all bound state poles deriving from the S-matrix; -ffei,...,en 
are normalisation constants, and Xi = 

The only unknowns are therefore the functions (5ei,...,e„(a;i, . . . , x„). They are ho- 
mogenous functions, analytic apart from the origin and polynomially bounded because 
of ( |1.6| ). They are symmetric functions in their arguments Xj and Xj if = e^. 

Let us focus our attention on the fundamental particle of the theory. In sections |^-^ 
we will deal with systems containing only one asymptotic state and therefore only form 
factors with all indices corresponding to this particle have to be considered. Denote the 
corresponding function Qii...! = Qm for n particles. This is, according to the above dis- 
cussion, a totally symmetric function. Therefore we introduce the elementary symmetric 
polynomials, which form a base in the space of symmetric polynomials and are generated 
by i 

n n 

Ylix + Xi) = J2x''-''akixi,X2,...,Xn) . (1.12) 

1=1 k=0 
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Further we will use the notation that x = 1/x and . . . , z„) = (Ti(xi, . . . Then 

the functions Qn can be cast into the following form 

wherein 7 denotes a possible non-integer power which is determined from the spin of the 



corresponding operator through (|1.4| ) and ( |1.5| ). The polynomials P„ can be expressed 
as polynomials in the elementary symmetric polynomials cxj. 

Before turning to the explicit construction of solutions of the form factor equations for 
specific systems, we want to discuss the general structure of the space of local operators. 
Suppose we have calculated the form factors of a scalar operator (for example the trace 
of the energy momentum tensor). Representing the operator in terms of its form factors 
one finds that the form factors of the derivative operators dzO ( dgO) are given by o'lJ-'n 
{d'lJ-'n)- This indicates, that in the massive model the two chiral sectors are represented 
by the variables Xi and Xi respectively. But since 

ai{xu...,Xn) = -, ^ , (1.14) 

any x dependence can be rewritten in terms of x, reflecting the fact that right and left 
chiral sectors are not independent in the massive theory. 

We want to introduce a grading into the space of operators that resembles the struc- 
ture of highest weight representations of the Virasoro algebra, in order to find their 
massive counter-parts. It can be seen from ( |1 . 1 1|) , that increasing the power in ( |1.13D 
for fixed spin s, will increase the divergence of the corresponding form factor JF„. For 
that we define as the chiral (left) operators those with the mildest ultraviolet behaviour, 
which means we put = in (|1.13 ). Similarly we conjecture that augmenting N in- 
creases the presence of right components in the derivative state. A simple example of 
this fact is = d(f) as we discussed above. Unfortunately we cannot assign a one-to- 
one correspondence between conformal states {h, h) and massive states (A^ + s, A^). The 
parameter A^ rather indicates how many operators L act on the primary operator, but 
also depends on the explicit representation of L„ in the massive theory. We will discuss 
this fact in detail for specific examples. 
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Finally we want to stress that setting the scale N is a rather subtle issue. Because in 
every module of chiral massive descendent operators, there will be polynomials P„ which 
are proportional to (7„. Therefore one could identify those as chiral descendents of the 
primary operator, and regard the original ones as corresponding to = 1 states. For 
some operators this ambiguity can be solved, since one knows the form factors of the 
primary operators explicitely. But as we will see in section ^ this scale is not necessarily 
the same for all primary operators in a theory. 

In the next section we analyse the scaling field theory near the critical point of the 
Ising model, in zero magnetic field, A^3,4 + 0i,3. Because of its simple structure (it can 
be described as a massive free fermion theory), we can introduce our method and test it 
against known results. In section |^ we investigate the Yang-Lee model, as the simplest 
example of an interacting field theory. In section ^ we analyse another non-unitary one 
particle model, A^3,5 + 0i,3. Finally, in section || we generalise these results to a series 
of models, namely A^2,2n+3 + 0i,3- As we will see the physical picture persists, while the 
mathematical methods become more involved. 



As a first step we will apply the form factor bootstrap to the scaling field theory near 
the critical point of the two dimensional Ising model in zero magnetic field. This massive 
field theory can be described in terms of free Majorana fermions, or equally in terms of a 
bosonic particle which interacts through an S'-matrix, S = —1. The theory is divided into 
two sectors. One contains the monomials in the fermion fields, and correlation functions 
of such operators can be calculated by using Wick's theorem. The other sector contains 
fields which aquire a phase e*'^ when moved around the fermions. The basic fields in this 
sector are the order and disorder fields and fi{x). Their correlation functions have 
been analysed in numerous ways, and using the form factor approach in @, 0]- 

For the scaling Ising model the parametrisation ( |1.11|) takes the explicit form 0, 0, ^ 



2 The Ising Model 



^„(A.../3„) = g 

n 5 • • • ; •^n )Hn{crn) 



1 



Y[ tanh 

i<j 



n 



(2.1) 




2 
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This parametrisation reduces the equation of kinematical poles 

Qn+2( X, Xi, . . . , Xfi) X Qni^-^ly ■ ■ ■ y •^n) ; 
Qn+2{~X, X, Xi, . . . , Xn) = , 1 



71) to 



1 
2 





(2.2) 



Notice the half-integer power of cr„ which appears in the sector w = 0. It has to be 
introduced in order to satisfy (|1.4|), and has the effect that form factors with odd index 
in this sector will carry half integer spin. 

The sector w = ^ contains the order and the disorder fields a and /i. The form factors 
of both operators have been calculated and are given by Qn = const, in ( p.l|) . In the 
descendent operators of cr{x) have been determined. They were constructed by building 



symmetric polynomials of spin s, which are invariant under the recursion relation (|2.2| ). 
Only invariants with odd spin are consistent with the recursion relation, and can be 
expressed in terms of the elementary symmetric polynomials as the n + 1 dimensional 
determinant [|I^ 

cTi (Ts 0-5 0-7 

1 (72 (T4 (76 
1 (72 
1 ^2 



'2n.+l 



(2.3) 



Arbitrary polynomials formed from the quantities ( |2.3| ) constitute solutions to 
which is the equivalent situation as in conformal field theory where the Virasoro irre- 
ducible representation of the operator 0i^2 is generated by the application of the operators 



-(2n+l) 



to the primary field. 



The method we want to describe chooses a different basis in the single subspaces of 
spin s. Namely we use the fact that the number of solutions of the recursion relations at 
level n is given by the number of solutions at level n — 2 plus the dimension of the kernel 
of the recursion relations ( p.2|) . This amounts in a simple counting procedure which only 
involves the comparison of the degree of the polynomials Qn and that of the kernel^. This 
has the advantage that one finds all solutions to the form factor equations, and not only 



^Note that the second line in (2.2) is precisely the definition of the kernel of the recursion relation. 
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those which can be considered as descendent operators of some non-trivial primary field. 
In the case of the Ising model the kernel of the recursion relation (|2.2|) is given by 



JCn = lliXi + Xj) . (2.4) 

i<j 

The total degree of this function is deg{lCn) = \n{n — 1). 

Comparing the degree of the polynomial Q„ and the degree of /C„, we can now con- 
struct explicitly solutions of the recursion relations. For the odd form factors, we find 
the following solutions for the first few values of spin s: 

s = Qi = const. , 

' = ' • (2.5) 

s = 2 Qi = al , 

s = 3 = af, Qs = JCs , 

which indicates that we have 1,1,1,2 linear independent solutions at spin levels s = 
0,1,2,3 respectively. Notice that these values correspond exactly to the first terms of the 
character expansion of the conformal field 0i,2, namely Xi.2 = 1 + q + q"^ + 2g^ .... 
We introduce the character x of the conformal minimal model as 

^M = qi^~^/^^)^M , (2.6) 

where c and h denote the central charge and the conformal dimension of the primary op- 
erator respectively. Therefore x encodes the degeneracies in the Virasoro representation 



11 



-| oo 

Further, define V{m, n) to be the number of partitions of n into numbers whose value do 
not exceed m. These partitions are generated by 

^ =^P(m,n)g" , (2.8) 



{q)m 

where 

m 

(g)^^n(l-?0 • (2-9) 



2=1 
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Using the definition (|2.8| ), we rewrite our counting of tlie solutions in (|2.5| ) in tlie 
following way. At s = the number of solutions is given by V{1, 0), for s = 1 by V{1, 1), 
for s = 2 by V{1,2) and for s = 3 by V{1,3) + P(3,0). Using this notation and the 
generating function ( p.9|) , one immediately finds the number of solutions at general spin 
s, as being generated by 

im(m,— 1) 

Fo^E ^j-^ = iSfiM) . (2.10) 
Note that this formalism leads rather to 'fermionic' expressions of the character than to 



the 'bosonic' form (2]7). This notation ( 'fermionic' and 'bosonic' ) has been introduced 
by R. Kedem et. al. [^] in their study of character identities for a large class of conformal 
field theories. Only recently bases have been constructed for conformal field theories 
which lead to fermionic sum expressions for the characters [|T^. It is an interesting fact 
that the form factor approach leads to them in a straightforward way. 

Similarly one can analyse the even form factors, being related to the disorder field /i. 
The counting procedure can be carried out analogously, and one finds that the number 
of chiral operators for spin s is generated by 

„im(r7i— 1) 

Go- E '-^r^ = xtt\<i) . (2.11) 

m,even \'djm 

As it is well known, both order and disorder fields are related to the conformal operator 

01,2. 

As a next step we investigate solutions corresponding to = 1 in ( |2.1[ ). We con- 
jecture that these operators correspond in the ultraviolet limit to states of the type 
Z_iL_„^ . . . 101^2) (compare figure ||). Let us repeat the counting as we did for the 
primary operators. The total degree of the polynomials Qn in this case is given by 
deg{Qn) = n + s which must be compared with the degree of the kernel of the recursion 
relations. 

We list the first few solutions, 

s = —1 Qi = const. , 
s = = ai, Qs = }C-i , 



s = 1 Qi = erf, Qs = aiJC 



(2.12) 



3 



s = 2 = af, Qs = crfiCa, Q3 = 0-2/03 
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Figure 1: Space of descendent operators of the irreducible Virasoro representation of the 

field 0i_2 of the conformal Ising model. 



The counting here gives not only the derivative states of the type L_iL„j . . . L„j,|(/)i^2) 
as expected, but also the original chiral ones. For example from ( |2.1j ) one immediately 
realises that the = 1 spin zero solution Qi = ai coincides with the = spin 
zero solution Qi = const., which is the one generating the form factors of the primary 
operator. In the same way every solution Qn being proportional to cr„ corresponds to an 
element of the set of chiral descendents. 

The generating function of the = 1 states can be found in an analogous way as 
before and is given by 

oo im(r7i-3) 
m,odd \H)m 

In order to find the level one descendent operators we need to subtract the chiral ones. 



rn,odd 



E 

m. odd 



oo „^m{m~3) f -\ /y™^ °° „4m(m— 3) 

{l)m m.odd il)m-l 



m,odd 



oo ^i(m+l){m-2) 

E — 7-. = qGo = QXiM) 

n,even \Q)m 
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where we introduced the notation q = q~^. We have carried out the calculation explic- 
itly, since this re-summation technique is a fundamental tool in the analysis of g-sum 
expressions and will be repeatedly used in the following. We find that in the perturbed 
model we have the same number of h = 1 + states as in conformal field theory. The 
analysis of the = 1 states in the even sector is in complete analogy, and will therefore 
be omitted. 

Finally let us discuss operators in the subspace N = 2, again for the odd sector. 
Generalising from the analysis of the subspace N = 1 one might expect that in the 
ultraviolet limit the operators for N = 2 correspond to states with h = 2 + j^. But for 
spin s = —2 we find two solutions, which we denote by 

gji) = const. , ^ = ai/Ca , (2.13) 

whereas from conformal field theory we would expect for h = 2 + j^, s = —2 only one 
operator (compare figure |l]). We conclude that the simple interpretation of N h does 
not hold in general. The further we increase N, the more terms the specific subspace will 
contain {e.g. for = 3, s = —3 we find 6 operators). This fact can be understood from 
the explicit representation of the invariants (|2.3| ). They have partial degree dp{l2n+i) = 
n + 1, or specifically for our example from above dpll^) = 2. This means we can write 
/3/al = Is, where /„ denotes the same invariants ( |2.3| ) but in the variables x. Using these 
expressions the states (|2.13| ) can be related to conformal states as 



Similar for all spins s in the subspace = 2 we expect states of the kind L^^L^m ■ ■ ■ -^-n^ 
This is indeed the case, since the generating function for the N = 2 solutions is given by 

im(m,— 1)— 2m 

= E —T-^ = {l + q + f + <f)XiM) , 

m,odd ^'i)">- 

that is, it resembles exactly the structure of Xi,2(^)Xi,2(q') up to second order in g, but 
additionally contains one term (fxi,2{l)- 

The presence of such higher level operators makes a direct link between the parameter 
N and the dimensions of the conformal fields difficult. Nevertheless increasing we can 

12 



gradually build up the whole space of operators in the massive theory. The number of 
operators at level is generated by 

oo ^m{m—l)—Nm 



m=l,odd 



in the odd sector, and by 



oo g^m{m—l)—Nm 

Gn= E 73 ' 

m=0,even \H/m 

in the even sector. These functions satisfy the recursion relations 

Fn = Fn-i — q^Gn-i , Gn = Gj^^i + q^FN_i , (2-14) 

with the initial conditions Gq = Fq = Xi,2('?)- Because of the symmetry of the recursion 
relations and the initial conditions, it follows that F^ = G^. This implies that F^ = 
(1 + g^)F/v-i which can be solved to give F^ = Gn = n^i(l + q'')Fo, which resembles 
the character expansion of Xi,2{q)Fo up to order A^. Note that the whole space of states 
is generated by limTv^oo Fn- This can be written as 

^oo = XiMxiAl) , (2-15) 

and similarly for Goo- 

A few comments follow, in order to clarify the meaning of these expressions. We find 
that Foo and Goo generate the space of operators in this sector. They have the identical 
form of the conformal expression of the space of descendent operators of the primary 
field 01^2- Further note that in (p. 15 ) the characters of the two chiralities are separated. 



This is a main feature of conformal field theories, where the right and left movers can 
be treated independently. Here in the massive theory the two chiralities are present, but 
q and q are not independent variables, and therefore such a factorisation can only be 
regarded as a formal expression. 

Finally we want to stress the fact that one obtains for finite some 'finitized' expres- 
sions for the Virasoro characters. In this approach they appear naturally in the counting 
procedure of states. Note that the same expressions turn up in the study of corner trans- 
fer matrices where the finitization is related to the fact that one has a discrete system of 
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finite size (see e.g. |]T6|)- It would be interesting to understand whether there is a deeper 
reason for this relation. 

The sector w = Q has been discussed in [|T7|. In this sector the kinematical recursion 
relation maps the form factors onto zero, J^„(/3 + in, (3, Pi, . . . ,(3n) = 0. It follows that 
form factors with different particle number n are not linked, and therefore any kernel 
solution represents an acceptable form factor from the point of view of the form factor 
equations. The analysis can be carried out in a similar way as for the sector w = \. The 
result is that the even form factors generate the descendent spaces of the energy operator, 
which scales in the ultraviolet limit as (|, |) and of the identity operator (0, 0). The odd 
form factors generate the spaces of the fermions (|, 0) and (0, |). The ultraviolet limits of 
the correlation functions can be calculated explicitly for these operators and the scaling 
dimensions coincide with the corresponding conformal operators. 

The thermal Ising model has been examined also by various other means. But, since 
it can be described as a free field theory, many methods and results remain confined 
to this special model. The interesting fact about the method we have presented in this 
chapter is that it can be straightforwardly applied also to interacting theories. In fact 
any massive integrable field theory having a description in terms of a scattering matrix 
can be investigated in this way. We will give several examples in the following sections. 



3 Space of Operators in the massive Yang-Lee model 



The scaling Yang-Lee model can be described as a massive integrable field theory [|T8 



Its on-mass-shell spectrum consists of a single massive scalar particle of mass m. Also 
the conformal field theory A^2,5 is quite simple containing only one primary field /i = — | 
besides the identity operator. 

The scattering matrix of the massive particles is determined by the following two- 
particle amplitude |ll8[ 

tanh ^(p — Id) 

with a = Ivr [|. 



■''We keep whenever possible the notation general since we can use then some of the results we will 
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The form factors have been analysed in |jTO|, |T9|, By defining the minimal two 
particle form factor as 



B { dr coshrf^ - ^1 sin^ ^('^-/^) 



(3.2) 



the parametrisation ( |1.11|) of the form factors reads for this model as 

(3.3) 



where Qn is of the form (|1.13| ). The constants if„ in ( p. 3D are given by 



a 



Hn = Co{4:Cos — sina {({i'K)^) 



2 



(3.4) 



The form factor equations have been reduced through this parametrisation to two 
coupled recursive relations for Qn [|^, namely 



Qn+2i, -^1 -^1 -^li ■ ■ ■ 1 -^n) Dfi{x^ • • • ; Xn)Qni^X\, • • • ; ^n) 

with the recursion function 

n k 



(3.5) 



(3.6) 



k=l m=l,odd 

wherein the symbol [/] is defined by 



[I] 



sin(/a) 



sm a 

The second recursion relation derives from the bound state axiom (|1.9| ) and reads as 



(3.7) 



Qn+l{xUjKxUJ 2,Xi, . . . ,X„_i) =X YI{X + Xi)Qn{x,Xi, . . . ,Xn-l) 



(3.8) 



i=l 



The recursion relation (|3.5| ) has been studied in |^ for the primary massive operators. 
It admits an infinite number of solutions. A basis in the solution space is given by the 
polynomials 

Qn = Qnik) = \\M,,{k)\\ , kez , (3.9) 

obtain here for the analysis of the form factors of the models Al2,2n+3, which will be discussed in section 
||. This is due to the fact that the S'-matrix of the fundamental particle of these theories is given by 
^]) with a = 2^'^- 



15 



where Mij{k) is an (n — 1) x (n — 1) matrix whose entries are 

Mij{k) = a2i-j [i - J + k] . (3.10) 

Enforcing in addition to (|3.5|) also the bound state recursion relation ( p.8| ) only one 
consistent solution exists, which is given by Qn(l) |0- In the ultraviolet limit the 
operator corresponding to this solution scales as the conformal primary operator 0i 2 
of the minimal model A^2,5 0]- ~ i the equation ( |3.5| ) admits only one single 



solution which is not consistent with the bound state equation 

We are interested in the analysis of the space of descendent operators of this model. 
From ( |1.13| ) we find that the functions Qn take the form . . . , x„) = (j~^P„(xi, . . . , x„) 



We start by analysing the space of chiral descendents, = 0. Note, that in this model 
even and odd form factors are coupled due to the bound state equation ( |3.8|) . The method 
is in principle the same as in the Ising model. We solve the relations ( p.5| ) and ( |3.8| ) re- 
cursively using the property that the space of solutions at level n is given by the number 
of solutions at level n — 1 plus the dimension of the kernel of the combined recursion 
relations. In table |^ we have written down the degree for the polynomials P„ required 
from Lorentz invariance and the dimension of the kernel of the recursion relations. 

Let us now count the independent solutions for various spin levels. For s = 1 the only 
solution is that generated by Pi = o"i, which corresponds exactly to the level 1 descendent 
of the primary field of the Yang-Lee model. For spin 2 we have two possible solutions, 
generated by the polynomials 

Pl=Crl , P2 = /C2 . 

Similarly higher spin values can be investigated. From table ^ we find that a kernel 
solution for n = 3 occurs only for spin s > 6. Let us examine the solution space for the 
particular spin value s = 6. The solutions are 

Pi = at , 

P2 = Cr^/C2, 0"2CT^/C2, Crl)C2 , 

P3 = /C3 , (3.11) 
that is we have 5 independent operators. 
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We compare these values with the dimensions of the spaces of descendent operators 
in conformal field theory. The number of descendent operators is given by the character 
expansions, 



n 


rfe^(Q(f)) = in(n - 1) + s 


deglC = |n(n — l)n > 2 


1 


+ s 





2 


l + s 


3 


3 


3 + s 


9 


4 


6 + s 


18 


5 


10 + s 


30 


6 


15 + s 


45 


7 


21 + s 


63 



Table 1: Comparing the total degrees of the polynomials P„ for the case N = with the 
kernel of the recursion relations for spin s and level n. 



Xi,i = l + g' + g^ + g^ + 2g^ + 2g^ + 0(g^) , 

Xi,3 = l + g + g' + g^ + 2g^ + 2g5 + 3g6 + 0(g^) . (3.12) 

Summing up these values we find that the Yang-Lee CFT contains 1 descendent operator 
at spin 1, 2 operators at spin 2 and 5 operators at spin 6. We compare these values with 
the ones obtained in the perturbed model and find that they coincide. 

In table ^ we have collected the number of independent solutions of the recursion 
relations using the notation of partitions introduced in ( |2.8| ). The counting can now be 
generalised to arbitrary values of s. A new kernel solution will contribute, whenever 
s = n{n — 1). Therefore the number of solutions is generated by 

n=0 V^Jn 

where we introduce the general definition, 

oo n^—Nn 

Fn^T.^--— . (3.14) 

n=0 
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s = 




S = I 




s = 2 


P(l,2) +P(2,0) 


s = 3 


P(l,3)+P(2,l) 


s = 6 


P(l,6) +P(2,4) +P(3,0) 



Table 2: Number of independent solutions of the recursion equations for several spin 
values s. 



For the Yang-Lee model the Rogers- Ramanuj an identities (see e.g. [12, 13|) relate the 
so-called p8|] 'fermionic' sum expressions to the 'bosonic' form of [0, [1^ as 

1 



i^o(l-g2+5i)(l-g3+5^) 

oo 1 
~(2,5). . 1 



iS?'fe) = n7r— r 



F-1 



(3.15) 



By a re-summation one can show that 

Fi = Fo + F„i = xi,i(g) + xi,2(g) , 

which means that the number of chiral descendent operators at each level is the same 
in the conformal and the perturbed Yang-Lee model. It is interesting to note that we 
have found the sum of the two characters, and not the two Virasoro irreducible repre- 
sentations modules separately. This is expected from the physical point of view, since 
in the perturbed model a mixing of the two spaces occurs. This has been analysed by 
Zamolodchikov [|T|, by determining the conservation laws of the perturbed theory. The 
simplest such conservation law, is given by 

where T denotes the holomorphic component of the energy momentum tensor, which 
is an element of the Virasoro irreducible representation of the identity operator in the 
conformal model. This relation expresses the continuity equation for the stress energy 
tensor in the massive model. 
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Next we investigate the solutions corresponding to Q„ = -^Pn, i-G. = 1 in (|1.13| ). 
Let us repeat the counting as we have done for the primary operators. The total degree 
of the polynomials P„ is now given by deg{Pn) = — 1) +n + s. The first few solutions 
are 

s = — 1 Pi = const. , 

(3.16) 

s = l P, = alP2 = ai/C2 , 

s = 2 P^ = af,P2 = af}C2, P2 = ■ 
As in the Ising model this amounts in an over-counting of the states, since the chiral 
operators appear in the subspace iV = 1 as solutions proportional to i-e. for which 

Pn ~ 0-n- 

The generating function of the = 1 states is given by F2. In order to find the 
structure of this subspace, we subtract the chiral operators using the g-sum relation 



F2-Fi= qFo = qxn{q) 



The structure of the A^ = 1 space coincides with that of the descendent operators of the 
conformal scalar primary fields and 0i^3 for h = 1. In CFT there are only states 
proportional to Xi,3{(l) since = 0. We find therefore that the subspace A^ = 1 in 

the massive model corresponds to the space h = 1 + Jiq of the conformal model, where 
ho stands for the right chiral dimension of the conformal primary fields. 

Using the generating function ( |3.14| ) one can show that the number of operators in 
the subspace A^, is generated by Fat+i. This function satisfies the recursion identity 



Fn — Fn-i + ^Fn-2 



(3.17) 



with the initial condition of Fq = Xhsil) and F_i = Xi,i(q'). It has been shown by I. 
Schur (see e.g. [O) that the solution to this problem is given by 



F^ = Df-'\q)xiA<l) + D'2'-''m 



1,1 



(3.18) 



with 



00 

E 



A^+ 1 



(3.19) 
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N+1 



1 + [ ^+1-5m ] 



(3.20) 



Z}W(g) = ^ (_l)A'^t(5/.-3) 

which are finitized forms of the characters Xi,3(^) ^ind Xi,i(^) respectively. From the 
recursion relation ( |3.17|) it follows that Fn = FN^i{modq^), and therefore we find that 
F/v resembles the partition function of the Yang-Lee model up to order g^. Similarly as 
in the Ising model we can write the formal identity 

hm F(iV) = ^2,5 = xMxiAq) + XiMXiAq) ■ (3.21) 

There are several new features in this model with respect to the Ising model. First 
in the form factor approach we cannot disentangle members of the Virasoro irreducible 
representations of the identity operator and the primary field. As we have explained this 
is due to the fact, that a mixing of these Virasoro irreducible representations occurs in 
the massive theory. On the other hand if the structure of the spaces are isomorphic one 
should be able to find a unique mapping of the operators. 

Actually a similar structure appears in the Ising model in the sector w = 0, which 
contains the descendent operators of the massive scalar primary fields 0i,3 and (pi^i. Also 
there the same argument of the mixing of operators holds. The special fact in the Ising 
model is that for the sector w = the UV dimensions of the operators can be calculated 
explicitly. Therefore it is possible to assign in a unique manner the form factors to specific 
conformal operators. It has been shown in [jl^ that in the Ising model {w = 0) chiral 
operators in the descendent space of the identity operator have form factors which are 
proportional to cr„ while the others form the chiral descendents of the primary field 0i^3. 

It is tempting to conjecture a similar structure in the Yang-Lee model. This is due 
to the fact that the space of solutions with P„ 7^ an are given by 

qn{n-l) 

E73 — = xi,i{q) ■ 

n [Dn-l 

This seems to be the opposite situation as in the Ising model, but one has to recall that in 
the Yang-Lee model the field 0i,3 has negative scaling dimension, and therefore behaves 
less singular as the identity operator in the UV limit. Unfortunately we are not able to 
calculate the UV dimensions exactly in this model, and therefore this conjecture cannot 
be confirmed. 
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Finally we want to discuss the relation of the space of operators to the conserved 
charges of the model. In [|1^ P. Christe has constructed descendent operators of the 
primary field. His approach was to modify the construction of ^ for the Ising model 
to the case of the Yang-Lee model, i.e. he constructed polynomial invariants which 
satisfied both recursion relations (|3.5| ) and (|3.8|) . He found invariants of spin values 
s = 1, 5 mod (6). They correspond to the infinite set of integrals of motions of the massive 
theory. Note that in the Ising model polynomials formed the corresponding invariants can 
be mapped one-to-one to the fields of the conformal Virasoro irreducible representation. 
Obviously this is not the case for the Yang-Lee model, since the operators of the Virasoro 
irreducible representation of 0i^2 are generated by L^n+i, L^n+i- Constructing derivative 
operators through polynomial invariants one only finds a sub-space of the descendents of 
the conformal operator 01^2- This implies that the conserved charges do not generate the 
whole space of operators. This feature generalises to all further models we will consider. 
The Ising model is an exception since it can be considered as a free field theory. 



4 The model 5 

As a further simple model whose on-shell spectrum contains only one particle is the 0i 3 
perturbation of ^^3,5. Let us recall the basic properties of the conformal model. It is a 
non-unitary minimal model, with the primary operators given through the Kac table 



3 

4 





1 


1 


5 


20 


1 


1 


20 


5 





3 
4 



(4.1) 



The characters of these operators are |]TT], P^ ], 

"("+2) 00 
n=0,even v"/" ^=—00 

71=0, even y^J'n- /x=— 00 
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n — I — 

n=0,odd WJn ^=-00 

n(n+2)-3 ^ 

XM(g) = E = E 9'''''-'^ - q^'^^'^^"-^'^ . (4.2) 

n=0,odd /i=-oo 

We have collected the bosonic and the fermionic expressions for later convenience. 

Perturbing the conformal model in the direction of the 0i^3 operator one obtains a 
massive field theory, containing one asymptotic particle without a bound state. The 



interaction of the asymptotic states is given through the S'-matrix |£6 



1 / TT 

S{P) = —i tanh 2 ~ ^2 

An interesting novel fact with respect to the previous models is, that these particles have 
non trivial statistics, since lim^^-too S{l3) = =Fz, and are in fact related to the operator 
02,1 = 01,4 of the model p^. Because they are spin j particles they also have Wij = t\ 
depending which chirality of the operator we choose (or equivalently if we define Wij 
through ±00 in rapidity space). 



The minimal two particle form factor can be determined as p5| 



^'"-(Z?) = sinh ^ TT -) ^ ) ^ . (4.3) 



Since there are no bound states the parameterisation of the form factor is given by 

^„(A, ...,/?„) = K^Qnixu ...,x^) a.Vn^'™'" n } ' (4-4) 

,-y» _l_ rf> 

where we define the normalisation constants Hn such that they satisfy 

The fractional power of cr„ in ( |4.4| ) in the sector = is needed in order to satisfy (|1.4| ). 
As a consequence we will find that operators corresponding to odd form factors will have 
non-zero spin, namely the same as the asymptotic particles, s = 

Using the parametrisation ( [4.4| ) and taking into account the statistics of the asymp- 
totic states, the recursion relation ( |1.7D reduces to 

Qn+2{-X, X, Xi, . . . , X2n) = (-x)*^""'" -D„(x, Xi , . . . , X„) Q„ (Xi , . . . , X2n) , (4.6) 
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where 



. . . = 5](e("-^)- - e2^™e-("-^)-)x"-V^ . (4.7) 

k=0 

We first analyse the sector w = 0. Even and odd form factors are not hnked by the 
recursive equations, and we start by examining the even ones. As before we use the 



parametrisation ( |1.13|) to introduce a grading in the space of operators. We have 7 = 



since the fractional powers of cr„ have already been taken in account in ([4.41). 



Since we have carried out the counting technique already in detail for the Ising and 
the Yang-Lee model we will merely state the results. For chiral operators (A^ = 0) the 

2 

degree of P„ is given by deg{Pn) = \ + s. For spin s = there exists only one solution 
which was found in p5|, and corresponds to the perturbation of the scalar primary field 
01^3. Its form factors are given by Qn = | |cr4j_2j-i| |, n = 2,4,... in ([4. 41), where the 
dimension of this determinant is |. 

For spin s = 0, 1, 2 one finds 1, 1, 3 solutions respectively, while for general spin values 
the generating function for the number of independent solutions is 

n(n-2) 

n,even \Hjn 

By shifting the summation index and using the fermionic sum expressions of the charac- 
ters ('14.21) it is not difficult to show that 



i^r = xi,i(g) + xi,3(g) . (4.9) 

We want to analyse the full space of operators, and therefore increase the value N in 
( |1.13|) by integer powers. The generating function for states at level is given by 



n(n+2) 

4 



-Nn 



Let us further define 



n,even vi)'' 



Nn 

(7 4 "'^ 

g%^j: ^-j--- , (4.11) 

n,odd vyjn 



then these functions satisfy the recursion identities 

^N - -^N-l + Q '^N-l 



G% — G%_i + Ff 



(4.12) 

N 
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In order to separate formally the chiral sectors, we express these infinite sums in terms 
of the characters. Prom the recursion relations it is clear that and can be written 
as polynomials in q multiplied by the characters, which are the initial conditions of the 

recursion relations, namely Fq = Xi,i{q) = Xi,3{q)- Expanding for the first few 

values of N one finds indeed that 



The general solution of the recursion relations can be cast into the form 



(4.13) 



with 

D'^\q) 



/i= — CX) 



D'^\q) 
D^'\q) 



lx=—oo 



-(3/i+l)(5At+l) 



_ -(3,.+l)(5;i+2) 



-M(15/i-4) 



/*=— oo 



■(3m+1)(5m+3) 



lx=—oo 



r/i(15//+7//+l) 



_ ^(3m+2)(5m+1)+1 



2N 
N -l-bfx 

2N 

2N 
N -2-bii 

2N 



N -2-5n 



XiM > (4 

^ xi,3(g) , 

^ q.xiM) 



2N 
N -5fi 

2N 
2N 

2N 
N -l-bji 
where the limits are understood as — * oo. 

While Ff^ generates the even part of the partition function as we expect, it seems 
that the function G% lacks a physical interpretation in this context. Its importance will 
be clarified in the following. 

Let us introduce a compact notation in order to identify the operators corresponding 
to a specific set of form factors. If the space of operators in a certain sector can be 
decomposed as Xr,s{Q)Xr',s'{Q) we will call the corresponding chiral {i.e. with the lowest 
possible divergence of the form factors) operator with spin s — hj. g — hj.f gf , as the per- 
turbed primary field {hr,s,hr',s')- In this notation the space of even form factors in the 
sector w — contains the descendents of the primary fields (|, |) and (0, 0). 
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The analysis of the even form factors is complete, and we turn to the odd sector. 
From ( |1.5|) and (|4.4| ) we find that it will contain particles of spin — ^ + Z. By applying 
the counting method we find that the space of chiral descendents is given by 

n(n-2) 

E = Q~hi,2{q) + qhiAiq) (4.15) 

n,odd WJn 

which indicates that 0i^4 does not mix with the primary operator 0i 2 but with its de- 
scendent at level one. Compare this situation with the table of conformal dimensions 
of the model The dimension of the field 0i,4 {h = |) is by far 'closer' to that 

of L_i0i^2 {h = than to that of the primary operator itself {h = — ^). Therefore 
this result is not surprising. On the other hand this indicates that when studying more 
complex systems, where the dimensions of the primary operators differ by several units 
(as it is the case for example for the unitary minimal models), this counting procedure 
might become rather complex, since different primary operators will mix with descendent 
operators at various levels. 

For general we use again generating functions Fn and Gn now defined as 



Fn^E ^-4-. and . (4.16) 

n,odd y^Jn n,even \H)n 

Observe that they have the same functional form as the generating functions for the even 
sector (|4.10| ) and ( |4.11| ), but the summation over odd and even indices is interchanged. 



4 



They satisfy the same recursion relations ([4.12|) but with the initial conditions 



3 1 

q^XiA^q) ^"^^ ^0 — 9 *Xi,2- This implies that the space of operators can be decomposed 
as 

F^ = q'hi,2{q)xiAq) + qhiAiq)xiAq) ■ (4.i7) 

We have calculated the form factors of the primary operators explicitly, and they are 
given by the polynomials Qn in ( |4.4| ) 

Qn = \W4i-2j-i\\^~^ , n=l,3,..., (4.18) 

for the field (~2l)' |) where the upper index of the determinant indicates the dimension 
of the corresponding matrix. For the operator (|,0) the polynomials are 

n — 3 

Qi = ai and Qn = <yn\\o-4i-2j\\^~^ , n = 3,5,... . (4.19) 
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We now turn our attention to the sector w = ^. The recursion function ([4.7|) takes 
the form 

Dn = T. ^kx""-' COS ^{n-k) . (4.20) 
The counting argument determines the space of solutions, as being generated by 

^0= E ff = Xi,2(g) . (4.21) 

In order to understand the structure of the space of operators we also analyse explicitly 
the terms with = 1. It is given by 

'^1= E = + + XiAid) , (4-22) 

n,even 

that is 01,4 does not mix with the primary operator 0i,2 but again with its first descendent. 
We have drawn a picture of the perturbed space of operators in figure ^. 




Figure 2: Space of descendent operators in the massive model A^3,5 + 0i,3 for the sector 
w = \. 

The full space of operators is generated by which becomes in the limit n ^ oo 

Glo = xi,2(g)xi,2(g) + ggxi,4(g)xi,4(^) ■ (4.23) 
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The form factors of the primary operator (~^) — ^) can be calculated explicitly, and are 
given by 

Qn = \WAi-2j\\^^~^^ , n = 2,4, 



in ( [4 .41 ). In this example it becomes clear that the expressions describing the content of 



operators in a specific sector in terms of the characters have to be regarded as formal 
ones. Especially in ( |4.23| ) there appears the term qq which multiplies the characters of 
the field 0i.4. If we cancel them, we would not alter (4.23), but the expression would no 



longer reflect the actual structure of the space of operators in this sector. 

Finally let us summerise the analysis of the odd form factors in the sector w = ^. We 
find that the space is structured as 
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Gl, = q-^XiMxiM + Q *Xi,i(9)Xi,4(g) • (4.24) 
The form factors of the primary fields are 

Qn = \\o'4i-2j-l\\ ^ , ?T,= 1,3,..., 

for the operator (|, — ^) and 



^ n — 3 \ 

Qi = ai and Qn = \\o-4i-2j+i\\ ~ , n = 3, 5 



for the operator (0, |). 

Let us discuss the physical implications of the above analysis. The massive model 
7^3 5 + 013 contains the following primary operators: the scalar fields (|, |) and (0,0) 
and the para-fermionic fields ( s = — ^ + Z), namely (— |), (|; 0) in the sector w = 0. 
In the sector w = | we find the fields (~^)~^)) (f;!) and the spin s = ^ + Z operators 
(|, — ^), (0, |). We want to stress the fact, that we are not able to establish a direct link 
between the form factors of these massive operators and their conformal counterparts, 
but have rather defined those operators through the degeneracies of their descendent 
spaces. 

An interesting feature of this model is the appearance of the para-fermionic operators. 
Their spin is determined by the spin of the asymptotic states. This causes that only 
specific combinations of chiral dimensions can appear in the massive model. Thus in 
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contrast to the space of operators in the conformal theory where arbitrary combinations 
of right and left dimensions can occur for a theory defined on the real line. 

We have introduced the chiralities in the massive model through the variables x 
and X = K Therefore, in order to have a consistent description we should be able to 
transform the form factors of the operator (~2l)' |) ^'^ those of the operator (|, — ^) by 
the transformation x —>■ x (and similarly for the other pair of para-fermionic fields). Let 
us analyse the form factors of (—^5 1), given by 

^^^■^^ = ff„a|-^||a4.-2,||<"^> n f l i ■ (4-25) 

yXi ~r Xj J 



We define J-'n{Pi, • • • , Pn) = ^n{—Pi, • • • , —Pn)- Using the relation (|1.14|) we find for the 
determinant 



n-l 



|cr4i-2j||^"2^^ = a„ 2 ||cT4i-2i~i||^"2'> , (4.26) 



and therefore 



jr( 20.5) _ jr^5' 20) _ (4^27) 



- (0 -) (-0) 

Analogously one can show that JF„ = jFn*' . 

Related to this fact is the following observation. In calculating the form factors, we 
had to choose which chirality i.e. which limit of the S'-matrix to select in the analytic 
continuation process. We chose the limit P ^ 00. Suppose we would choose S{—oo) = +i. 
Repeating the same analysis as before, one finds that the scalar fields remain in the 
respective sectors, while the spin-dependent fields change place. This is as expected, 
since now the chirality of the spin is now linked to the right component of the field 
dimensions. 

A final remark concerns the scalar operators of the theory. We have found (pi^i and 
01^3 in the sector w = while we found 0i,2 and 0i,4 in the sector w = It is interesting 
to note that these locality properties are in accordance with the operator-product expan- 
sions of the correspondent conformal operators with the field 0i 4. It has been observed 



previously that the field 02, 1 (which for ^^3,5 is the same as 0i,4) is deeply connected 
to the operators generating the asymptotic states. This fact is not totally understood. 
The analysis of more complex theories in the formalism proposed in this article might 
give new insights to this problem. 
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5 The models M2 2^+3 

From our point of view, the models 2,2^+3 constitute the simplest systems containing 
more than one asymptotic particle, since their Kac table contains only one row of fields. 
Also their spectrum contains only scalar particles, and from the analysis of the previous 
sections we expect only scalar operators to appear in the massive model, as it was also 
the case for the Yang-Lee model. 
The S'-matrices are given by |]23[ 



min{a,b) — l 

Sab = f\a^b\^{P)f(a+b)^{P) U (/(|a-M+2fc)f (Z^))' , (5-1) 

k=l 

where a = and a,b = 1,2,. . .p label the particles of mass = sin (af)- The 

functions / are given by 

tanh ~ ^Q^j 

Note that the S'-matrix of the fundamental particle has the same functional structure as 
the S'-matrix of the Yang-Lee model, the difference being that the parameter a is specified 
differently. We can therefore use the results from section |^, where we had analysed the 
kinematical recursion relation for the form-factors corresponding to this S'-matrix element 
for general a. 

The form factor equations for the models A^2,2p+3 + 4>i.3 have been analysed in detail 
in There the form factors for all perturbed primary operators (pi^k, k = l...,p 

have been calculated explicitly. They are given by ( ^TSf ) with Qn = Qn{k) where Qn{k) 
is defined in (|3.9|) . Further it was shown that there is only one independent recursion 
relation resulting from all possible bound state fusion processes. This recursion relation 
can be expressed as a constraint on the form factors with only indices 1 ... 1 relating 
F„+2p+i to Fn. It follows that the consistent form factors of this model form a subset of the 
solutions of the kinematical recursion relation on particle 1. This allows us to determine 
the structure of the form factors of higher particles by carrying out the bootstrap on the 
expression obtained for -Fi i i. 

We carry out the bootstrap explicitly for the form factors involving particles one and 
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two. They are given by 



-^2,2,...,21,l,...,l(^i'---'^9+-) = (5-3) 



m=l i<j m=l i=l m<l 

The single components in this complex expression are 

^ (^.)sinh(f )C„W _ 

[Xi + Xj) < a > 

" te„i.C(!ir,a)2 < 1 - f X fa > ' '^'^^ 

4xmX„7r('(z7r, a)"' < 2a >< a >< 1 — a > [Xm + Xn) 
Hn is the normalisation constant of and 

A{x) = — cot q; tan— — . (5.7) 

Y n({iTr, 2a) C(^7r, a) 2tt 

Further we have introduced an abbreviated expression for the functions introducing the 
bound state poles, 

< a >= sinh -{Pij — ia) sinh -{Pij + ia) , (5.8) 

and the minimal 2 particle form factors 

Cii(/5) = C(^,a) , 

Ci2(/5) = C(^^,f)C(^,f) , (5.9) 
C22W) = aO,2a)aO,aY , 

with 

f. r{k + ^ + i;-f)r{k + i-i;-f)r(k + l + i; + f)r(k + 2-i; + f 
L{ r(fc + i-|L_|!)r(fc + iL_M) r(*+|-t + ¥)r(^ + i + |, + f 

We start our investigation by analysing the system Ai2j- The above expressions are suf- 
ficient to examine the operator content of this model, since it contains only two particles. 
Let us analyse the form factors corresponding to chiral operators, i.e. = in ( |1.13| ). 
This implies that the functions Q2,2.. .2,1,1. ..i are polynomials. 
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dq = 



do = 12 



c'q = 3 \ dQ = 5 

nil \ H211 \ H221 \ ^r222 

dye = 9 dye -11 d^c = 16 I dye - 24 



Figure 3: Comparing dimensions of Qei,...,e„ with the respective kernels of the recursion 
relations. The lines indicate the inter-relations due to the bound state equation ( tL.9|) . 

As before, we apply the counting technique, but now we have to take in account that 
we treat a two-particle system. The recursion problem for the functions Qei,...e„ is now 
of matrix- form. We need to determine the respective degrees of the polynomials Qei,...en 
and of the kernel of the recursion relation. The degrees of Q2,2...2,i,i...i follow from Lorentz 
invariance, as 



deg{Q2, 2 2,1, 1 i) = -(2g + r)(2g + r - 1) - g + s 



(5.10) 



where s denotes the spin of the operator as usual. The degrees of the kernels of the 
recursion relations, /C2,2,...,2,i,i,...,i follow from the pole structure of ( |5.4j) - (|5.6| ), and are 
given by 

degi)C2^ 2, . . . , 2,1, 1, . . . , l) = - 1) + 4rg + 4g(g - 1) . (5.11) 

q r 

We have depicted these dimensions in figure |^. As in the one particle systems we 
use this comparison of degrees in order to count the number of independent solutions of 
the recursion relations. Note that the polynomials Qei,...tn symmetric in coinciding 
indices only. Therefore in this example we have two sets of indices, which have to be 
taken in account separately. From figure ^ we immediately find that the first solutions 
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are being generated by 

W)i W)i 

which correspond to the possible choices of one-particle form factors. The terms 

+ + 7^ , (5.13) 



{q)2 {q)i{q)i {q)2 

generate the number of independent solutions for the two-particle form factors. The power 
of q in the nominator derives from the spin- value at which deg{Qei,e2) + s = deg{)Cei^e2)- 
In particular, the denominator {q)i{q)i corresponds to Q2i{xi, X2) and indicates that this 
function is not symmetric with respect to its two arguments. Finally for the three-particle 
form factors we find 

+T^ + T^ + ^ ■ (5-14) 



(9)3 {q)2{q)i {q)i{q)2 (g)3 

In order to extract the general generating function, one uses the fact, that the spin 
value at which a new solution enters is given by deg{ICej^^,,,e„) — deg{Qe^^,,,^^). Defining 

„(ni+n2)^+n|— a(ni+n2)-bn2 

Fa,^T.- T-rj-. , (5.15) 

ni,n2 \H)rii\H)n2 



and using (|5.1CI|) and ( ^.11|) one finds that the space of chiral operators in the 0i^3 per- 



turbation of M.2,7 is generated by Fi^i. The functions ( |5.15| ) satisfy 

Fa,b = Fa,b-i + (f^'^ , 

= + g"-iF„_2,b . (5.16) 

These recursion relations can be used to prove that 

^1,1 = ^-1,-1 + ^0,-1 + ^0,0 = Xi,i(g) + Xi,2(g) + Xi,3(g) ■ (5.17) 

As in the previous cases we find that the chiral operator content of the massive model is 
in a one-to-one correspondence with that of the conformal model. 

Unfortunately the results for this model are less complete than the others. We have 
tried to solve the recursion relations (^.16|) to find an exact expression for the whole space 



of operators (given by \\mN~*oo Fn,n) of the model, but did not succeed. Nevertheless we 
have solved them numerically up to order 50, and the results indicate that 

Fn,n = xi,i(g)xi,i(g) + xi,2(g)xi,2(g) + xiAq)xi,3iq) ijnodq^) , (5.I8) 
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as expected. 

Finally let us turn to the general systems A42,2p+3- Again we need to determine the 
dimensions of Qei,...,en and the kernel /Cei,,,,,e„. Denote as Q"^i.™-2,. -,mp ^j^g function Qei,...,en 
corresponding to rrii indices e = i, where i = 1, . . . ,p, the number of different asymptotic 
states in the theory. From Lorentz invariance one finds that 

^ i=i 1=1 1=1 ^ 

For example the dimension of Quip^ii'^i)-, i-e. the function corresponding to the two- 
particle form factor of particles k and /, has dimension deg{Qki) = kl . 

In order to determine the dimension of the kernel we must carry out the bootstrap 
on the form factors of particle 1, satisfying the kinematical recursion relation. This is 
straightforward for the denominators and for the factor sinh(|) in ( |5.4| ). But one needs 
to be careful, since further zeros appear from the bootstrap of the factor (u. For example 
for the two particle form factor J^ki with k > I one finds through the bootstrap that 

Ci,i- (sinh^^ "\(/3,(fc-0f) c(/5,(fc + /)|) nc(/5,(fc-/ + 2z)|) X 

''-^^-^ /a \ 

nnc /5,(fc+0iT-(^+j> ■ (5.19) 

i=2j=0 V ^ / 

The first line in ( p.l9| ) gives exactly the minimal two-particle form factor for the particles 
k and l,Cki, while the second line reduces to trigonometric functions, due to the relation 
({P, (y)C{P, —d) =< >. They cancel poles which arise from the bootstrap of the poles 
in the denominator of ( |5.4| ). 

We will not give a detailed account of this calculation. Carrying out the bootstrap on 



all factors in (|3.3|) one finds that the total degree of the kernel of the recursion relations 
corresponding to the function Q^i is given by degiJCki) = k.l + 2.1. Comparing with the 
dimensions of Qku one finds that solutions corresponding to the kernel of K,ki will enter 
at level 21 (recall we have the convention that k > I). Therefore the general exponent in 
the g-sum expression will be given by 

2lnkni for k>l , 
knk{nk — 1) for k = I , 
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for every pair of indices Uk and ni. Taking the sum over all particles we find 

Nk , (5.20) 



k=l \ l<k J k=l 



where we defined Nj = X^iLo '^p-* • This leads to the final expression for chiral massive 
states as 

N^-Ni+N^-Ni+.-.+N^-Np 
mt^np {<i)nA<i)n2 ■■■{<!) np 

We want to show that this sum expression coincides with a sum over characters of 
the model. For that we need the g-sum expressions of the characters, given through the 
Gordon-Andrews identities (see e.g. [|^) as 



Xi,,(9) = i^00... 0-1-1. ..-1 > (5-22) 



where we have introduced the generating functions 

q 

which satisfy the recursion identities 



Nl-kiNi+Nl-k2N2 + ...+Nl-kpNp 

Fki,...,kp^ E ^:w3 > (5-23) 



Fki,...,kp -ffci,...,fci_i,fci-l,fei_|-i+l,fci+2v,fcp~'~ 

q'^-'q'^-' . . . q'-^-'Fki-2,k2-2,...M-2M+i,...,kp • (5.24) 
Apply these identities to the generating function of the chiral operators, namely 

-^11...! = -^111. ..10 + -^-1-1. ..-1-1 

The last factor corresponds to Xi,i(^)- we have to continue with 

-^11. ..110 = -^11. ..01 + -^-l-l,...-10 

This last term does not correspond to any conformal character, but if we use the recursion 
relations to permute through the to stand at the first position we find Xi,2{q) along with 
a plethora of terms, without physical interpretation. If one carries out this commutation 
process of generating always further zeros as the last index and commuting them through, 
one will finally end up with 

i^ii...i = Exl?"^H^) ■ (5.25) 



1=1 
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We conclude that in this series of models the operator content of the massive theory 
coincides with the content of conformal scalar operators. Unfortunately in all these 
systems, containing more than one asymptotic state, we were not able to solve the general 
recursion- relations ( p. 24 ) which determine the full space of operators of the massive 
model. We have done several numerical calculations which lead us to the conjecture that 
the full space of operators in the massive model corresponds to the scalar conformal one, 
i.e. 

Z = i2xi,imiM ■ (5-26) 

i=l 

We have also compared the structure of the solutions of these recursion relation with 



several known 'finitized' expressions for the characters of these models [^, but were not 
able to find one which provides a solution to the recursion problem ( p. 241 ). 



6 Conclusions 

The aim of this work was to analyse the structure of the space of operators in massive 
integrable systems. Our particular interest was in deformed minimal conformal field 
theories, because these theories at the critical point have a well-known structure, that 
is the operators are organised into families determined through the Virasoro irreducible 
representations, with a highest weight state corresponding to a primary operator with 
rational conformal weight. All other operators can be obtained by applying the generators 
of the Virasoro algebra to these primary fields. 

For massive integrable models such an algebraic analysis has not yet been developed. 
As an alternative we have used the form factor bootstrap method in order to determine the 
structure of the space of operators. In order to find similarities we have used the physical 
properties of the conformal operators and have searched for analogues in the massive 
models. In particular we have introduced a grading into our space of operators in the 
massive model defining chiral operators as those with the mildest ultraviolet behaviour 
and through their spin. The whole space of operators is constructed by augmenting the 
divergence of the form factors gradually. 

Using the recursive equations for the form factors we have proposed a counting proce- 



35 



dure in order to find the number of operators in eacli of the finite dimensional subspaces 
defined through this grading. Though we have analysed only simple systems (in the sense 
that the Kac-table of the corresponding conformal model consisted at the most of two 
rows), the method is applicable to any system described through a factorised scattering 
theory, including systems with degenerate particles and even massless scattering theories, 
which describe the flow from one critical point to another. 

Our analysis of the massive models gives several interesting results. For all the sys- 
tems examined the chiral field content is in a one-to-one correspondence with that of the 
underlying conformal theory. The massive primary operators are grouped into families 
according to the parity of the corresponding conformal field. This induces the possibility 
of a mixing of the descendent operators in the corresponding Virasoro irreducible rep- 
resentations in the process of perturbation. The full space of space of operators in the 
massive model contains two types of fields. For theories with scalar asymptotic states 
there are only scalar fields, while for models with asymptotic particles which satisfy 
non-trivial statistics, in addition to the scalar fields also para-fermionic operators are 
present in the theory. This generates for example in the model A^a.s operators satisfying 
para-statistics with fractional spin |. 

Finally let us comment on the mathematical methods involved. A feature, which 
persisted in all models examined, was that the counting of the operators lead naturally 
to fermionic sum expressions of the characters, or rather sums of them. Recently a lot of 
interest has evolved in this subject, after Kedem et.al. [EH] have conjectured a whole series 



of such character identities, which have been proved in [g^, by using finitizations 

of the character expressions. As we have seen, the counting procedure leads to finite 
expressions for the characters, when we generate the full space of states. It seems to be 
an interesting fact, that the method here proposed gives a natural physical application 
to these mathematical results. 
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